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ABSTRACT 

Aims. It is known that stellar activity may complicate the analysis of high-precision radial-velocity spectroscopic data when looking 
for exoplanets signatures. We aim at quantifying the impact of stellar spots on stars with various spectral types and rotational velocities 
and comparing the simulations with data obtained with the HARPS spectrograph. 

Methods. We have developed detailed simulations of stellar spots and estimated their effects on a number of observables commonly 
used in the analysis of radial-velocity data when looking for extrasolar planets, such as radial-velocity curves, cross-correlation func- 
tions, bisector velocity spans and photometric curves. Stellar and spot properties are taken into account, as well as the characteristics 
of the spectrograph used (generally HARPS). The computed stellar spectra are then analyzed in the same way as when searching for 
exoplanets. 

Results. 1) A first grid of simulation results (radial- velocity amplitudes, bisector velocity-span amplitudes and shapes, and photom- 
etry) is built for F-K type stars, with different stellar and spot properties. 2) It is shown quantitatively that star spots with typical 
sizes of 1% can mimic both radial- velocity curves and the bisector behavior of short-period giant planets around G-K type stars with 
a vsin; lower than the spectrograph resolution. For stars with intermediate vsin/, smaller spots may produce similar features. Such 
spots may complicate the search for low-mass planets on short-period orbits. In these cases, additional observables {e.g., photometry, 
spectroscopic diagnostics) are mandatory to confirm the presence of short-period planets. We discuss these possibilities and show 
that, in some cases, photometric variations may not be enough to clearly rule out spots as explanations of the observed radial-velocity 
variations. This is particularly important when searching for super-Earth planets. 3) It is also stressed that quantitative values obtained 
for radial-velocity and bisector velocity-span amplitudes depend strongly on the detailed star properties, on the spectrograph used, 
on the line or set of lines used, and on the way they are measured. High-resolution spectrographs will help in distinguishing between 
spots and planets. 
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1. Introduction 

Most of the exoplanets detected so far have been identified 
thanks to the so-called radial velocity (RV) technique, and many 
more are expected to be found with presently ongoing surveys. 
Moreover, thanks to precision improvement, planets with masses 
as low as a few Earth masses are a priori detectable, and it is ex- 
pected that planets with masses similar to the Earth's will be 
detectable in the future if precision of a few cms ' is reached. 
However, it has been known for a long time that stellar activity 
(spots, pulsations) can also induce RV variations that can be pe- 
riodic and mimic those induced by planets. Then, there might be 
risks of misinterpretating RV variations when they have periods 
shorter than or equal to the star rotational period. 

An illustration of potential difficulties encountered with 
stellar phenomena is 51 Peg i tself. |Gray| (199 7) argued that 
the RV variations reported by Mayor & Queloz| (1995) and at- 
tributed to an orbiting planet with a 4-day period were instead 
due to some stellar phenomenon. [Gray] (1997) indeed claimed 
the presence of bisector shapes characteristic of stellar spots 
or pulsations. However, no quantitative modeling of the am- 
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plitude expected for RV or bisector variations were provided. 
It rapidly turned out that these claims were wrong: in fact 
[Hatzes, Cochran and Bakker| (1998) did not confirm those pe- 
culiar bisector variations on the basis of more data at a higher 
signal-to-noise ratio (SN). Also, very precise photometry did 
not reveal any variations down to a level of A(b + y)/2 = 0.2 ± 
0.2mmag ( Henry et al.| 1997), sufficient enough to reveal the 
presence of stellar spots. Finally, the planet around 5 1 Peg was 
confirmed, even though lBrown et aL] (1998), based on a detailed 
analysis of the possible effect of pulsations, concluded that they 
cannot exclude the possibility of pulsations as the source of 
51 Peg RV variations. This example shows that stellar activity 
(spots, inhomogeneities, pulsations) has to be taken into account 
when analyzing RV data of potentially short -period planets. We 
focus in the following on star spots. 



When modeling stellar activity, ISaar & Donahuel (1997) 
made first quantitative estimations of the impact of stellar spots 
on the RV curve of Fe I lines at about 6000 A, in the simple 
case of an equatorial spot with Tspot = K on an edge-on solar- 
type star They showed that peak-to-peak RV amplitudes up to 
a few hundred ms"' can be produced by spots or convective 
inhomogeneities, depending on the spot size and the projected 
rotational velocity of the star. They derived quantitative laws for 
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the RV amplitude, as well as for the bisector velocity-sparj]] am- 
plitude. They showed that the effect of the star's projected rota- 
tional velocity is important for detecting RV variations, but even 
more for detecting bisector variations. |Saar, Butler and Marcy| 
(1998) also estimated the impact of inhomogeneous convection 
from the study of the bisector velocity-span itself. They conclude 
that inhomogeneous convection also leads to RV and bisector 
velocity-span variations up to a few tens ms"^ They measured 
the weighted dispersions on RV measurements on a sample of 
G and F-type stars and found that, for G-type stars with ages 
about 0.3 Gyr (with rotational velocities of about 8-10kms~'), 
cr[, ranges between 20 and 45ms"'. Later on, Paulson et al.l 
(2004) measured a jitter up to 50 ms"' due to stellar activity in 
a sample of Hyades dwarfs. Note that lSaar & Donahuel (1997) 
had also pointed out that convective inhomogeneities can lead to 
even larger RV variations, especially for G2V-type stars. 

Generally, observers have tried to examine either the pho- 
tometric curves or the bisector variations, in addition to RV 
data, whenever RV variations with short periods were observed. 
As far as photometry is concerned, the maximum amplitude of 
variations in milli-magnitude is ^ 2.5/ (Saar & Donahue 1997), 
where / (in %) is the fraction of the stellar disk covered with the 
spot. Depending on spot size and photometric precision, such 
variations may or may not be detectable. An example of the 
impact of these complete studies is HD 166435 (- 200 Myrs, 
GOV), for which Queloz et al. (2001) rejected one short-period 
(3.8 days) planet candidate, based on bisector measurements, Ca 
II lines, and photometric observations. 

In this paper, we investigate in more detail the impact of stel- 
lar spots on RV curves and on other diagnostics that are com- 
monly used to disentangle cases of stellar activity from those 
of planets. We compute the visible spectra of stars with vaii- 
ous spectral types (from F to K), projected rotational velocities 
and orientations, covered with one spot of different sizes and at 
different positions. Then, we quantify the resulting RV, bisector 
velocity-span, and photometric variations. We finally present the 
results and discuss the impact on RV studies. 



2. Description of the simulations 

We define the following parameters for the star: temperature Teff , 
gravity log^, rotational velocity Vrot, inclination /, metallicity 
[Fe/H], microturbulent velocity Vfnicro macroturbulent velocity 
Vmaci-o- We assume Vmicro = 1.5 km s"' and Vmacro = 0.9 km s"' for 
a G2V star, and a limb-darkening coefficient e = 0.6. For the 
spot, the free parameters are the spot colatitude on the star 
surface, the temperature Tspot, and the spot size described by the 
parameter fr, which is the fraction of the visible hemisphere cov- 
ered by the spoj^ 

2.1. Spectrum computation 

Our simulations use Kurucz models (IKuruczl l993). The 3D stel- 
lar surface is divided into longitudinal and latitudinal sections. 
The number of cells is tuned so that the velocity sampling is bet- 



ter than the resolution of the generated Kurucz spectrum. The 
resolution of Kurucz spectra is chosen to be twice the intrinsic 
resolution of the instrument we wish to simulate. 

For a given set of stellar and spot properties, we first com- 
pute a synthetic spectrum using Kurucz models without rota- 
tional broadening. Then, we apply this spectrum to each cell, 
shifting the spectrum according to the Doppler law, taking its 
radial velocity into account. We then sum up the contributions 
of each cell, weighted by the cell's projected surface and limb- 
darkening, and taking the presence of a spot into account, when 
relevant. If a cell is covered by a spot, we use the black-body law 
for each wavelength of the spectrum to evaluate its weight com- 
pared to the same cell without spot {i.e., at Tspot = T^-ff). We gen- 
erate the resulting stellar spectra at different epochs of the star's 
rotational phase. Typically, we take 20 epochs to cover the phase. 
Each spectrum is convolved with the instrumental point spread 
function (PSF). The spectra are computed in the range 377- 
69 1 nm, corresponding to the High-Accuracy Radial-velocity 
Planet Searcher (HARPS, Mayor et al.| 2003) wavelength range 



' The bisector velocity-span measures the global slope of the bisector 
( IHatzesl l996) and allows shape variations of the lines to be detected. 

^ For a small spot (or <c I) defined by its semi-angle a (2a is the 
angle under which the spot is seen from the star center), /j- = I - cos a. 
Note that this fraction is not identical to the fraction of the projected 
area covered by the spot fp on the 2D stellar disk, which is equal to 
sin^ a for a low spot. For small values of a, we have = 2/,-. Hereafter 
the spot size used is assumed to be /- unless specified. 



or, in some cases, in a range corresponding only to one order of 
the spectrograph (see below). 

2.2. RV, bisector, and photometric variations 

We select, for each order, ranges for the spectrum that will be 
considered to compute the RVs and cross-coiTelation functions 
(CCFs). This selection is made so that we include neither broad 
lines, such as those of H or Ca, nor those from telluric origin. 
From these synthetic spectra, we then used the method described 
in|Chelli|(2000) and lGalland etan (2005a) to derive the RVs, the 
CCFs, bisectors, as well as bisector velocity-spans for each order 
of the spectrum or for the global spectrum. The method for com- 
puting the RVs consists in a kind of correlation in Fouiier space 
of each spectrum and a reference spectrum built by summing 
up all the spectra, which are specific to this star This method is 
fairly equivalent to the CCF one for G-K type stars, but is more 
suited to the case of A and F-type stars with vsin/ typically 
greater than 10-15 kms"'. The bisector shape of the lines and 
the resulting bisector velocity-span are estimated on the CCF in 
the same way as in Queloz et al. (2001). Note that, whenever we 
wish to compare with real data, it is possible to add noise to the 
spectra so as to mimic real data as closely as possible. Finally, 
we also derive the photometiic curve over the rotational period 
at 550 nm. 

In the following we compute the RV, bisector, and photo- 
metric variations over a wavelength range naiTower than the full 
HARPS wavelength range (some orders were not even consid- 
ered in some cases), and it corresponds to the part of the spec- 
trum that is not affected by atmospheric lines. This choice was 
made in order to stick as closely as possible to the reduction of 
actual A-F data, presented for instance in lGalland et al.l (2005a.b 
and 2006a,b). 



2.3. Example 

Figure [T] shows an example of a spot located at a colatitude 9 
- 30°, covering 1.02% (/■) of the visible surface of a G2V-type 
star (T^ff = 5800 K), with vsin/ = 7kms"' and / = 30°. The 
RV curve, CCFs, and bisectors are shown, as well as the bi- 
sector velocity-span as a function of RV. Finally, we also show 
the photometric curve. In this example, we see that the spot 
produces an RV curve similar to what can be observed when 
a planet (0.51 Mjup with a 3.7-day period) orbits a star with a 
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Fig. 1. Example of a spot at 6* = 30° with a size of = 1.02% 
on a G2 V-type star rotating with v sin / = 7 km s"^ and seen with 
an inclination / of 30° with respect to the line of sight: star spot, 
RV curve, CCFs, bisectors, bisector velocity-span curve, photo- 
metric curve. Values measured on the whole spectrum. Part of 
the spectrum used: orders #10 to #58. A Keplerian model with 
0.51 Mjup and a 3.7-day period would fit the RV data well (bot- 
tom). 



peak-to-peak amplitude of 145 ms"'. The bisectors clearly pro- 
duce an umbrella-like shape, and the bisector velocity-span is 
clearly correlated to RV variations. Its peak-to-peak amplitude 
is ^ 120ms The photometric amplitude is 1% peak-to-peak. 
This case is similar to the one of HD 166435 ( |Queloz et al.| 
2001). 



The estimations of 5 show a higher dispersion; nevertheless, 
a general trend is also observed (see Table |A.l"] i. This higher dis- 
persion observed for the bisector velocity-span comes from the 
fact that CCFs and resulting bisector velocity-spans are more 
sensitive than the RVs to line shapes and blends, given our RV 
measurement method. As a matter of fact, the CCF is the aver- 
age of several lines with shapes that differ from one to another; 
averaging over a single order with a small number of lines (a few 
tens of lines per order presently) does not completely average out 
the effects of individual line shapes. This effect is even stronger 
when considering stars with larger v sin /, as more blends, even 
with faint lines, affect the line shapes. 

There are several consequences: 

- If we are satisfied with results on A with a 10% range of 
validity, then we can use only one single order If we need 
results with better accuracy, we have to compute the whole 
stellar spectrum. 

- Given the observed dispersion on the bisector velocity-span, 
quantitative results or laws derived from a single line or set 
of lines (e.g., [SaaF&Danahue, 1997 or lHatzesl 2002) should 
not be used for quantitative comparison with "global" CCFs 
(i.e., obtained using the whole spectrum). 

- Inter-comparison of S values based on averaged CCFs over 
the whole spectral range should be restricted to stars with 
similar spectral types and projected rotational velocities be- 
cause of their strong effect on the broadening of lines. 

- The chromatic dependence of A, even though limited to a 
level of about 10%, may be used as a very precious addi- 
tional criterion for spot /planet diagnostics. This will be de- 
veloped more later Also, observing in near-infrared will al- 
low reduction of the effect of spots and thus allow detection 
of less massive companions of active stars. 

Note that, even though these trends have been shown in the case 
of an edge-on solar-type star, they would qualitatively remain for 
other geometries. 



4. Edge-on G2V-type star - equatorial spot 

We have so far computed a limited series of simulations of com- 
plete spectra for different values of star's projected rotational 
velocities, inclinations, and different spot characteristics (see 
Table [T}. In this section we present the results for an edge-on 
star (/ - 90°) with an equatorial spot (0 - 90°). 



Table 1. Cases for which full spectra were used (G2V-type 
stars). 



3. Order-to-order variations 

We first tried to check the dependence of the obtained values 
of A (peak-to-peak amplitude of RV variations) and S (peak-to- 
peak amplitude of the bisector velocity-span variations) on the 
spectral orders. Figure |2] shows the values obtained for A and 
S in the case of vsin/ = 2, 3, 5, and 7kms"' as a function of 
the order (i.e., of the wavelength). We see that the estimations of 
A vary within about 10%, with a smooth trend, decreasing with 
increasing order. The slopes of the trends for the different v sin / 
are given in the Table A. 1 We see that the absolute value of the 
slope increases with v sin /. 
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Fig. 2. Measured RV amplitude (top) and bisector velocity-span 
amplitude (bottom) as a function of order number; assuming an 
equatorial spot with /,. - 1.07% size on a G2V-type star seen 
edge-on and with vsin/ = 2, 3, 5, and 7kms"' (from bottom to 
top curves). Part of the spectrum used: orders #10 to #58. 
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Fig. 3. A 1.07% equatorial spot on an edge-on G2V-type star 
rotating with vsin / - 2kms"': star spot, RV curve, CCFs, bi- 
sectors, bisector velocity-span curve, photometric curve. Values 
measured on the whole spectrum. Part of the spectrum used: or- 
ders #10 to #58. In this case with a vsin / value lower than the 
spectrograph resolution, the bisectors are only shifted (without 
change of shape) even if the origin of RV variations is the pres- 
ence of a spot. 



4.1. General results 

For a given spot location, the values of A, S , and the bisector 
shapes strongly depend on the star's projected rotational veloc- 
ity. Figures [3] and |4] show as an illustration the RV, bisectors, 
bisector velocity-span, and photometric curves obtained respec- 
tively for vsin / = 2 and 7kms"'. The peak-to-peak amplitude 
of the RV increases from 37 to 155 ms i.e., roughly propor- 
tional to vsin/. This was expected from the first-order simu- 
lations quoted above. More striking and important are the dif- 
ferences observed in the bisector variations. Even though we 
observe an umbrella-like shape for the bisectors in the case of 
high-enough v sin / (note that this umbrella-like shape is also ob- 
served in Fig.[T]), as well as a clear correlation between A and S 
(inclined "8" shape), we only observe a shift in the case of low 
V sin /. In this case, when the v sin / value becomes lower than 
the spectrograph resolution (non-resolved lines), the effect of a 
change of shape in the stellar lines only results in a shift (at first 
order) after convolution with the instrumental PSF, which one 
can understand easily: the shape then results essentially from 
the shape of the PSF. In such a case, it will not be possible to 
find correlations between the bisector velocity spans and RVs. 
We estimate that the astrometric effect of this 1% equatorial spot 
on an edge-on star at 550 nm is a photocenter shift of 73 fiAU 
(thus, 7.3 juas at 10 pc), assuming a stellar radius equals 1 R0. 



4.2. Impact of spot size and v sin / on the RV amplitude 

We first checked with a full spectrum that the effect of the spot 
size is the same for every order To do so, we compared the A val- 
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Fig. 4. Same as Fig. [3] but with vsin / = 7kms"'. In this case 
with a V sin / value greater than the spectrograph resolution, the 
bisectors undergo changes of shape from the presence of a spot 
that also leads to RV variations. 
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Fig. 5. A (filled circles) and S (open circles) as a function of the 
spot size /,., for an equatorial spot on an edge-on G2V-type star 
with y sin / = 7 km s"' . Values measured on the whole spectrum. 
Part of the spectrum used: orders #10 to #58. 
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Fig. 6. A and 5 as a function of the star's projected rotational ve- 
locity, for an equatorial spot on an edge-on G2V-type star Values 
measured on the whole spectrum. Part of the spectrum used: or- 
ders #10 to #58. A varies approximately linearly with vsin/. S 

varies as ^ -^(v sin 0^ + Vg - vq j . If we take into account the in- 
strumental resolution (dashed line) a ^ 1.8, otherwise (dotted 
line) a - 2.5 (see text). 



4.3. Impact of spot size and v sin / on the bisector 
veiocity-span 

Similarly, we checked that the S/fr values do not depend on 
the spot size. Then, using spectra over the whole order range, 
we find that in the investigated range of values of vsin/ (2- 
lOkms '), S strongly depends on the star's projected rotational 
velocity, as can be seen in Fig. |6] More precisely, we find that 



S = 6.5 fr ^ ^(v sin i)^ + - voj 



1.8 



(2) 



where vq is the instrumental width (vq fixed to 3kms"' given 
a resolution of 100000 in the case of HARPS). The resolution 
plays a crucial role for the bisector determination: if we do not 
take it into account, we obtain S - 0.79 /,. (v sin /)^-^, and the 

value for the fit is multiplied by 3. This means that the PSF 
affects the shape of the star spectrum, and the less the lines are 
resolved the stronger this effect. If v sin / is negligible in compar- 
ison to vo (lines not resolved), then S is negligible (in compari- 
son to A for example), which is consistent with the fact that the 
deformation of the stellar spectrum then only results in a small 
shift of the observed spectrum. 

Besides, we note that for a given v sin i,A/S does not depend 
on /, . Hence A/S depends almost only on the star parameters. 



spot 



and spot location. 



4.4. Impact of tlie instrumental resolution on the bisector 
velocity-span 

To test the impact of instrumental resolution, we computed sim- 
ulations of an equatorial spot (fr - 1%) on a G2V-type star, seen 
edge-on, taking the instrumental profile into account or noj^ 
When taking it into account, we assumed the resolution to be 
100000 (HARPS) or 50000. Two cases of stellar projected ro- 
tational velocities have been considered: 3kms"' and 7kms~^. 
Figure[7]shows the obtained values of A and S , order per order. It 
can be seen that taking the instrumental PSF into account leads 
to lower values of both A and S , and the values decrease when 
the spectral resolution decreases. We also see that the impact on 
S is much stronger than the impact on A. In the case of the bisec- 
tor velocity-span, the values obtained are a factor of two smaller 
than the ones obtained when neglecting the instrumental profile. 



ues obtained for each order, assuming a star with edge-on spots 
with sizes of 0.5 and 1.07%. The values are found to be identi- 
cal from one order to the next one within less than 1%. We then 
performed several simulations on order #31, assuming different 
values of f,- ranging between 0.5 and 2%. The amplitudes of RV 
variations measured on order #3 1 are found to be proportional to 
the spot size (Fig.|5]l. Altogether, this shows that the RV ampli- 
tude measured on the whole spectrum is proportional to the spot 
size. 

Figure [6]provides the amplitude of the RV variations. A, as a 
function of v sin /, measured on whole spectra, for a spot size of 
1%. We find that in the investigated range of vsin/, A depends 
almost linearly on the star's projected rotational velocity. More 
precisely, with/, fixed to 1.07%, we find A - 17 (v sin /)'^ Then 
we can conclude that 

A = 16/,(vsin/)i\ (1) 
where /,• is expressed in percent, and v sin / in km s"' . 



4.5. Comparison with previous results 

Assuming an equatorial spot with a temperature of OK, on an 
edge-on, solar-type star. ISaar & Donahuel (1997) find that the 
amplitude of the RV variations of the 6000 A Fe I line, expressed 
in ms"', follows the law A^ 6.5/°'^v sin / where Aj = A/2, f 
is the fractional projected spot size (previously defined fp), and 
V sin / the projected rotational velocity, expressed in km s ' . Note 
that they did not take any instrumental resolution into account. 
They also derive the following law for the amplitude of bisector 
velocity-span variations, in ms"': A^^span ~ 0.1 l/"''(v sin /)■'•' 

with As, span = 5/2. 

Using the Ca I 6439 A line, and also assuming an equatorial 
spot with Teff - Tspot = 1200 K on an edge-on solar-type star and 
a spectral resolution of 0.043 A (R = 150000). iHitiesl (2002) 
find As = (8.6y- 1.6)/° '' andAs.span = (22- 16V + 3.3V^)f° '^ 



^ Note that even without taking the instrumental profile into account, 
the spectrum is still sampled over 4096 pixels per order, resulting in an 
instrumental resolution of 300 000. 
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Fig. 7. A and 5 as a function of order number when taking the 
instrumental profile into account (crosses: R = 50000; open cir- 
cles: R = 100000) or not (filled circles, corresponding to R - 
300000). Two cases are shown: G2V-type star with vsini - 
3kms"' (2 upper panels) and with vsin/ = 7kms"' (2 lower 
panels). In both cases, an equatorial spot with 1% size is as- 
sumed. Values measured from order #10 to #58. 



where / seems to be the projected spot size (fp). He clai ms that 
his results are compatible with those of lSaar & Donahuel (1997). 



while assuming a very different temperature for the spot. This is 
in fact surprising as spot temperature certainly impacts the RV 
and bisector velocity-span amplitudes (see below). We find the 
laws ([T} and Q for an equatorial spot on an edge-on G2V-type 
star, simulated with HARPS. They are obtained using the whole 
spectrum between orders #10 and #58. We recall that fr = fp/2. 

One has to be careful when comparing the results of differ- 
ent simulations. Indeed, as illustrated before, 1) the choice of 
different lines or wavelength ranges to measure the different pa- 
rameters (already highlighted by Gray 1982), 2) the spot temper- 
atures, and 3) the instrumental resolution, all have a significant 
impact on the values obtained for the bisector velocity-span. We 
find a quasi linear dependence of A as a function of / (Fig.[5]l, as 
found by the other authors and believe that the differences found 
between the power dependence on / or v sin / for the amplitude 
(power 0.9 for lHatzesl 2002 and power 1 for the present work) is 
probably not significant. 

Considering an edge-on G2V-type star with v sin / = 7 km s"' 
and Teff - 5800 K, an equatorial spot with /,■ = 1.07% (fp = 
2. 14%), and r^pot = 100 K, for a single line at 6006 A, we find A 
= 200ms-' and 5 ^ 238 ms^'. ISaar & Donahue, (1997) find A 



180 ms-' and 5 =268ms-^ for T, 



spot 



■ K, and the results are 



in reasonable agreement. If we take another line at 6439 A with 



pot 



100 K, we find A = 192 ms-' and S = 125 ms-'. This 



shows that S strongly depends on the line chosen. If we now try 
to compare our results with those of lHatzes1 (2002) with the same 
parameters, except T^pot - 4600 K and a single line at 6439 A, 
we find A = 121ms-' and S - 81ms-', while he finds A = 
232 m s ' and S - 284 m s ' . The discrepancy is quite important 
there. It shows that probably the ISaar & Donahuel (1997) and 
IHatzes (2002) results were in fact not compatible. We note that 
using Hatzes (2002) laws with f - fr instead of fp would lead 
to more similar results (A = 125 m s-' and S - 152 m s-'). 

All this shows that one should not make any quantita- 
tive comparison of bisector velocity-span values or AjS ob- 
tained on real data with formulas such as those proposed by 
ISaar & Don'5iuel (1997). IHatzesl (2002). or this paper Dedicated 
simulations corresponding to the same instrument and star prop- 
erties should be performed, at least to calibrate the dependence 
for a given instrument. 

5. Spot at different latitudes on an edge-on 
G2V-type star 

We ran several single-order (#31) simulations for a spot of vari- 
ous sizes {fr — 0.5-2%), located at different colatitudes: 9 - 90, 
60, 30, and 10°. We find again that, for a given star, the AjS 
ratio does not depend on the spot size and that A//,- and 5//, 
follow laws as a function of v sin with powers similar to those 
obtained in the case of an edge-on star with a similar velocity. 

Globally, the amplitudes of A and S decrease with increasing 
latitudes of the spot. In those circumstances, more cases occur 
where bisector velocity-span variations may not be detectable. 
Of course, in the present case (edge-on star), the RV curve and 
the photometric curve will be constant half of the time, and it 
will be easier to distinguish RV variations due to a companion 
from RV variations due to a stellar spot. Nevertheless, cases with 
several spots might be problematic. 



6. Spot on an inclined G2V-type star 

Obviously this case is more interesting in the framework of RV 
planets as the spot may be seen over the whole rotational phase 
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Fig. 8. Spot located at 6* = 60° with a size of 1.03% on a G2V- 
type star seen almost pole-on (/ = 10°), rotating with vsin/ - 
2 km s"^: star spot, RV curve, CCFs, bisectors, bisector velocity- 
span curve, photometric curve. Values measured on the whole 
spectrum. Part of the spectrum used: orders #10 to #58. In this 
case with a low v sin / value, the bisectors are only shifted, the 
RV curve mimics a companion, and the photometric variation is 
small. A Keplerian model with 0.1 Mjup and a 4.4-day period on 
a circular orbit would fit the RV curve. 



of the star As before, we computed a few simulations over the 
whole spectrum (see Table [TJ, as well as more numerous single- 
order (#31) simulations of a spot on an inclined star. For the 
detailed simulations, we assumed a spot of 1%. 



6.7. Examples 

Figure [T](G2V-type star, / = 30°, = 30°, 1.02% spot, vsin/ = 
7 km s-^ and Fig.]8](G2V-type star, / = 10°, = 60°, 1.03% spot, 
vsin / - 2kms"') show examples of results obtained when the 
whole spectral range is used (orders #10 to #58). Table |2]gives 
examples of values obtained for A and S using the whole spectra 
and a spot with /,• ^1%, and the diagnostics that can be used to 
detect if it is activity-induced RV. A/S is independent of fr, thus 
those values can be scaled down. We note that the features pre- 
sented in Fig. |8] are identical to those expected from an orbiting 
planet (0. 1 Mjup with a 4.4-day period). Orbiting planets produce 
periodic RV curves. The CCF shape does not change with time; 
instead, the CCFs are just shifted one to the other, consistently 
with a Doppler-shift of the spectra induced by the presence of a 
planet. 

In fact, for stars with non-resolved lines (low v sin /), when 
relatively low-amplitude periodic RV variations with short pe- 
riods are observed, the bisectors are just shifted and similar to 
those produced in the case of a planet. Additional observables 
are mandatory for attributing those variations either to planets or 
to spots. 
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Fig. 9. RV-curve shapes for various configurations (see text). 



6.2. Discussion of various cases 

In the case of low v sin /, for a given star inclination, the de- 
tailed shape of the RV curve depends on the spot latitude. 
Figure |9] shows different shapes that can be obtained with vsin/ 
- 2kms"'. We can distinguish several cases: 

- A: star at high inclination (/ = 90°). The spot is hidden part 
of the time, its projected size varies considerably along the 
rotational phase (equal to when the spot is not visible), thus 
photometric variation is high (Table |2]i. 

- B: star fairly inclined (/ = 30°). The spot is seen along the 
whole rotational phase (with = 40°), and its projected size 
does not change as much as in the previous case. Photometric 
variation is smaller (Table [2|. 

- C: star very inclined (/ - 10°). The spot is not very far 
from the equator (0 = 60°), it is seen along the whole ro- 
tational phase, but its projected size variation is smaller than 
in the previous case. Photometric variation is much smaller 
(Table|2ll. 

- D: star very inclined (/ - 10°). The spot is close to the pole 
(0 - 10°), it is seen along the whole rotational phase, and 
its projected size varies very little along the rotational phase. 
Photometric variations are very small (Table [2]i. 

Apart from the first curve (A), we can find a Keplerian model 
that is acceptable, even if those with small / or with a spot near 
the pole are better. For intermediate inclinations, the departure 
from a 1 -planet Keplerian model is small and the Keplerian fits 
are quite acceptable (Fig.[T]l. 

Concerning the RV and bisector velocity-span amplitudes, 
it can be seen that, for a given v sin /, A and S decrease when 
the star and spot inclinations decrease. Moreover, A and S have 
different values depending on v sin /. More precisely, we can dis- 
tinguish several case^ 

- For large vsin / (around 7kms"'), A/S ~ 1.0-1.5; hence, 
when RV variations are detected, bisector velocity-span vari- 
ations will also be detected in most cases, but for very small 
spots where it is below the detection limit while A is still 
detected. At high inclination, we may find cases where vari- 
ations in A are detected, whereas those of S are not {e.g., /, 
= 0.1%, / = 90°). However in such cases, variations will not 



We assume that A or 5 will be detectable with good confidence if 
they are larger than 6 ms"' (with HARPS); we recall that A and 5 refer 
to peak-to-peak amplitudes 
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Table 2. Values obtained for A (ms"'), S (ms"') and AV 
(mmag) for a G2V-type star with vsin/ = 2, 5, and 7kms"', 
and diagnostics that can be used to detect if it is activity-induced 
RV (v: RV, in some cases; b: bisectors; p: photometry). 
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A/5 
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10 
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48/27 


76/75 
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P? 


b,p? 
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be observed over the whole phase, so there is no risk of con- 
fusion between planets and spots, except in very particular 
and unlikely configurations with several spots. At very low 
inclinations (/ ~ 10-20°) and low 6, A/S ~ 1, which means 
that, if RV variations are detected, variations in S with simi- 
lar amplitudes can be detected. 

- For low vsin/ (2kms"'), A/S » 5, and it is independent 
of f/, S is always smaller than 3 ms"', whereas A is larger 
than 6ms"' for fr = 1 %. Moreover, in the case of low incli- 
nations of star and spot, the photometric variations may be- 
come smaller than 5 mmag, i.e., difficult to detect routinely. 
In the example shown in Fig.|8](/ = 10°, - 60°), the ampli- 
tude of photometric variations is 5.7 mmag (peak to peak), 
twice as small a spot would produce twice as small a photo- 
metric variation. 

- For intermediate v sin ; (5 km s"'), A and S variations are al- 
ways detectable for = 1%. A simple scaling shows, how- 
ever, that for much lower values of /,• (e.g., 0.2%), RV vari- 
ations could be detectable while S would not. This would 
happen for low values of / and 0, and the amplitude of photo- 
metric variations would become undetectable (at a 0.4 mmag 
level for a 0.2% spot). It is thus mandatory to carry out photo- 
metric measurements below the 1 mmag level to detect such 
problematic cases. 

- There still remains a narrow range of situations where A 
can be detected, whereas S cannot. This would happen for 
i ~ 30°, ~ 0.02%. In that case A ~ 3 ms-' and 5 ~ 
2ms'. And the corresponding AV would be on the order of 
0.3 mmag, hence undetectable. 

7. Equatorial spot on other types of stars 

We computed various simulations for two other types of stars, 
namely a F6V and a IC2V. Here we present the results of those 
simulations and analyze them. 



7.1. F6V-typestar 

We assumed an edge-on F6V-type star with a 1% spot. The 
star's projected rotational velocity was assumed to be 2, 7, or 
20km s"'. We assumed as before a difference of 1200 K between 
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Fig. 10. Measured A and 5 as a function of order number; as- 
suming a 1% equatorial spot on an edge-on F6V-type star with 
vsin / = 2 (squares), 7 (circles), and 20kms"' (triangles). 



the star and the spot temperatures: T^ff = 6200 K, T^poi - 5000 K 



We also assumed v„ 



1 .5 km s 



Figure 10 shows the values of A and S obtained for differ- 
ent V sin /. We find the following relations between A, S , fr, and 
V sin /, with vq still fixed to 3 km s"' (see Fig. 14 1: 



A = 15.4/,- (vsin/)' 



S =7.1/, (^(vsin/)2 + v2-vo) 



1.5 



(3) 



(4) 



In Fig. 1 1 we show an example of comparing with real ob- 
servations and the present simulations: the RV of an F7V-type 
star with vsin / - 9.6 km s '|^ measured with HARPS and our 
software, are variable (left, top), but the profile of the CCFs 
changes as a function of time (left, bottom). Our simulations 
of the same kind of star show that these variations (right, top: 
phased) can result from the presence of a spot induced by stel- 
lar activity: we obtain the same behavior of the CCF and the 
same correlation between the corresponding bisector velocity- 
span and RV (right, bottom), with the same amplitude levels. 
The value of logR'^^ is rather high, namely -4.3, which effec- 
tively indicates a relatively high level of activity. We note that 
this is coherent with the high level of amplitude variations (> 
200 m s"' peak to peak). 

7.2. K2V-type star 

We also conducted full-order simulations in the case of a K2V- 
type star with Tgff = 4800 K, Tspot - 3600 K, and Vmacro = 



5 Estimated from the CCF. IGlebockil (2000) gives V sin i = 7 km s , 
the simulation uses 1 1 km s"' to reproduce the slope of the correlation. 
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Fig. 11. Example of comparison of an active F7V-type star 
(vsin/ ~ 9.6 km s"') between real observations (left) and the 
present simulations (right): we obtain similar RV variation am- 
plitudes (top), CCF behavior (center) and correlation between 
bisector velocity-span and RV (bottom). 

0.9 kms '. The star was assumed to be seen edge-on and the 
spot to be equatorial. Different values of vsin / were used: 2, 3, 
5, and 7 km s"' . Figure 



12 



shows the values of A and S obtained 
for different v sin /. We find the following relations between A, 
S , fr, and V sin with vo still fixed to 3 km s"' (see Fig. 14 1: 



A = 18.3 /.(vsin/) 



1.1 



10.0 



V sin 0^ + ' 



(5) 
(6) 



As an illustration, we tried to find whether we could ex- 
plain the RV and bisector velocity-span variations observed with 
HARPS in the case of a KlV-type star. Figure 13 shows such an 
example: the RV, measured with HARPS and our software, are 
variable (top) with an amplitude of 130ms"', but the profile of 
the CCF changes as a function of time (center and bottom), with 
a bisector velocity-span amplitude of 50ms"'. We find that a 
close-to-equatorial, 1 % sized spot located on a star rotating with 
V sin i = 5 km s"' could produce similar amplitudes. Note that the 
value of log R'^^. is -4.3 for this star, which is consistent with the 
presence of a relatively high level of activity. 

7.3. Comparison to G2V-type stars 

Table [3] summarizes the values of A, S , and AV obtained for 
these edge-on F6V and K2V-type stars with various v sin / and 
a 1% equatorial spot. It allows for comparison with G2V-type 
star (second column). Figure 15 shows the comparison between 



F6V and K2V spectral types with respect to G2V order-to-order 
A and S measurements. Globally, one sees that, for a given set 
of conditions (spot location and size, star inclination, and v sin /, 
Jetf , and Tspot), the amplitude of RV and bisector velocity-span 
variations increase when we go from F6V to K2V-type stars. For 
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Fig. 12. Measured A and 5 as a function of order number; as- 
suming a 1 % equatorial spot on an edge-on K2V-type star with 
vsin / = 2 (squares), 3 (circles), 5 (triangles), and 7 kms"' (dia- 
monds). 
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Fig. 13. Example of comparison of an active KlV-type star 
(v sin / ~ 4.8 km s"') with real observations (left) and the present 
simulations (right): we obtain similar RV variation amplitudes 
(top), CCF behavior (center), and correlation between the bisec- 
tor velocity-span and RV (bottom). 
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Fig. 14. A and 5 as a function of the star's projected rotational 
velocity, for an equatorial spot on edge-on F6V, G2V, and K2V- 
type stars. A varies approximately linearly with vsin/. S varies 

as ^ -^(v sin i)^ + Vg - vq j . If we take into account the instru- 
mental resolution (dashed lines) a ^ 1 .5 for the F6V-type star, 
1.7 for the K2V, otherwise (dotted lines) a ^ 2.4, and 2.5. Values 
measured on orders #10 to #58. 

Table 3. Comparison between F, G and K-type stars. 
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a given v sin /, this means that it will be more difficult to detect 
bisector velocity-span variations in the case of F6V-type stars 



180 
160 

gl20 



■■ ■■■ • ■ ■ ■ 



^ " ■ AAA a- ■ -A- AAAAAAAAAAAA 
®00- £OOCO©3 



AAAA . A-aAA- 



-Gcoo- 



50 
40 
: 30 



10 20 30 40 50 60 70 

order 



200 
150 
100 
CO 50 



. • . -'^ ■ 

■■■ -A . ■ ■ 

«i ■ ■ 



A 

A -*AAAA^ ■ - -A 

„ AA . . _ -A*- -AS " A A 

AA .AaM- ■ AA A *0 

? I 



10 20 30 40 50 60 70 
order 

Fig. 15. Comparison between F6V (squares), K2V (triangles), 
and G2V (circles)-type stars for two vsin/ (open symbols: 
2kms-', filled symbols: 7kms-'). Orders #10 to #72. 



than in the case of G2V-type stars and easier for K2V-type stars. 
However, one has to note that F6V-type stars rotate an average 
faster than G2V-type stars. The occurrence of cases where possi- 
ble confusions between planets and spots should thus decrease, 
whereas K2V-type stars rotate an average slower than G2V-type 
stars, so the occurrence of those cases should increase. 



8. Impact on the search for short-period RV planets 
and further studies 

In the previous sections, we have shown that spots can produce 
a variety of features (RV, bisector shapes, and variations), whose 
characteristics vary according to the spots and star characteris- 
tics. We have also seen that the precise values of RV and, to a 
greater extent, the values of the bisector velocity-span depend on 
the line or set of lines taken into account and on the spectrograph 
used. A quantitative comparison of the bisector velocity-span or 
bisector velocity-span to RV correlation should be made for the 
same line or set of lines, acquired with the same spectrograph 
(or at least with the same spectral resolution) and analyzed with 
the same software. Working on the full CCF allows individual 
effects to be averaged out and allows quantitative comparison 
to models, provided the models assume spectral types and pro- 
jected rotational velocities identical to those of the star under 
study and use the same spectral lines and same spectral resolu- 
tion. 

From our studies, spots with typical sizes (1%) at the surface 
of stars with high v sin / will be easily identified and may be char- 
acterized using the criteria presented above (bisectors, bisector 
velocity-spans). In the case of stars with low inclinations and 
spots near the pole, confusion may however arise when the RV 
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amplitudes are small (< 10 m s"'). The situation is more complex 
in the case of stars with low v sin /. 

For stars with low v sin /, spot features may, in some cases, 
mimic those produced by planets. This happens when 1) the ob- 
served RV variations have periods similar to the star rotational 
period, 2) variations are observed in the RV curves, and at the 
same time 3) bisectors do not change in shape, but are just shifted 
according to the RV changes. From the study above, this happens 
in the case of stars with intermediate vsin/ (typ. 5kms"') and 
a very small (0.1%) spot or in the case of stars with low vsin/ 
compared to the spectrograph resolution (typ. < 3kms ' for a 
resolution of 100000, < 6kms"' for a resolution of 50000), 
even with a spot that has a common (typ. < 1%) size. When 
the spot is large enough and its location favorable, photometry 
will in general tell whether the variations are due to spots or not. 
In the case of small spots and for particular inclinations of the 
star, the photometric precision may be out of reach. Of course, 
the occurrence of such cases will increase when searching for 
planets with lower and lower masses (super-Earths or Earths) 
for a given orbit, as they will produce smaller RV amplitudes 
and bisector shapes comparable to those of smaller spots, i.e., 
with small S . Note that indicators of activity such as logR'^^ 
may not be sensitive enough for low levels of activity (hence the 
low amplitude of RV variations). Indeed, for levels of activity as 
low as logR'^j^ ~ -4.8, the level of the amplitude for the RV 
variations is still a few ms ' jSantos etal. 2000, |Wright et al. 
2003). 

The chromatic dependence of A may in such cases help in 
distinguishing between planets and stellar spots. Indeed, in the 
case of a planet perturbation, no such chromatic dependence is 
expected, whereas in the case of spots, chromatic effects will 
occur. Figure 16 for example shows the RV curves as a func- 
tion of order number in the case of a G2V-type star (v sin / - 
2kms"') seen with / = 30° and a 1% spot located at 6* = 30°. 
This spot would produce a photometric variation of 1.3%. The 
peak-to-peak difference between the measurements of A in the 
first 4 orders and the last 4 orders is 3.8 ms~', hence could be 
detectable with good signal-to-noise data. This provides a pos- 
sible additional criterion for testing the origin of RV variations 
when bisector velocity-span and/or photometric criteria cannot 
be applied. It has the great advantage of being an observable 
present in the spectroscopic data themselves. 

Will it be possible to invert an observed data set (RV vari- 
ations, bisector velocity-span variations, possible photometric 
variations)? The unknown parameters are /, number of spots, 
temperature of spot, spot colatitude, and spot size. We already 
know V sin / increases with earlier spectral types. One spot then 
results in 5 free parameters and 2 spots result in 8 free parame- 
ters, so it will probably be very difficult to perform this inversion 
reliably. Nevertheless, our study shows that, with a given set of 
observations, it is possible to make the diagnostic of the presence 
of spots, provided we are able to measure the chromatic impact 
on A. 

This study has assumed a very simple case of a single spot. 
Of course, the reality is more complex: the star may be cov- 
ered by several spots with different temperatures, at different lat- 
itudes; they may also have inhomogeneities distributed in com- 
plex patterns; in addition the star may undergo complex patterns 
of pulsation. Simulations will be performed in the future to ex- 
plore a wider variety of cases and to test more realistic cases. 
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Fig. 16. RV curves as a function of order number The high- 
est amplitude coiTesponds to the lowest (blue) order and the 
lowest amplitude to the highest (red) order. Here AV = 1.3%. 
Hypotheses: G2V-type star, vsin/ = 2kms"', / = 30°, 6 = 30°, 
/.= !%. 
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Appendix A: Linear fits 
List of Objects 



'51 Peg' on page 1 
'51 Peg' on page 1 
'51 Peg' on page 1 
'HD 166435' on page 
'HD 166435' on page 
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Table A.l. Values a and b of the linear fits y = ax + b corre- 
sponding to the different plots. 
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